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I. INTRODUCTION 

Random motion of a classical particle in a potential that breaks spatial inversion sym- 
metry and fluctuates in time generally leads to unidirectional flow. On the molecular level, 
many important processes in biology and nanotechnology rely on this mechanism. Biolog- 
ical motors produce mechanical work from metabolic energy in order to effect intracellular 
transport, or self-propulsion of bacteria through rotatory flagellar motion.- Another pro- 
cess ubiquitous in any living organism is active transport of molecules or ions against a 
chemical potential gradient. Such molecular or ionic pumps generally consist of a specific 
channel across a cell membrane. The internal binding sites of the channel are linked to con- 
formational fluctuations, stimulated by metabolic energy or by light.-"- Similar processes 
are known for artificial nanopores, with potential applications in molecule or ion separat- 
ing devices.-"^ In the realm of quantum transport, a net electron drift under an applied ac 
driving signal can be generated by various mechanisms, which may become analogous to 
classical unidirectional transport when dissipation is included.-i^ 

Our goals in this paper are twofold. First, we set up a model that describes unidirectional 
transport along a finite open system, represented by a non-symmetric discrete chain with 
time-dependent driving. This model emphasizes (i) coupling of the two chain ends to reser- 
voirs and (ii) interaction effects between transported particles that occupy the chain. These 
combined features distinguish our study from most works on Brownian ratchets, reviewed 
ii>i2iii. Indeed, both interacting Brownian motors^^i""— and ion pumping mechanisms that 
involve multiple occupation of the associated channel structure^ have become important sub- 
jects of research in nonequilibrium statistical physics and biophysics. The model we examine 
allows us to study production of both mechanical and chemical work. In what follows we 
refer to a machine working against a mechanical force as motor and to that working against 
a chemical load as a "chemical pump" .- Accordingly, attention will be focused on the system 
performance in both modes of operation and its dependence on the driving characteristics 
and on fundamental input parameters such as chemical potentials of the reservoirs, load 
potential, and interaction strengths. Although our model is certainly far from describing 
realistic systems, we argue below that at a qualitative level, the effects of concentration and 
interaction evaluated here should have rather general validity. 

Our second goal is to provide a test of time-dependent density functional theory (TDFT)i2i"— 



when applied to far from equilibrium dynamics under time-dependent driving signals.— This 
method can be regarded as a version of dynamic mean-field theory, constructed such that 
it can account for exact static properties. Its accuracy in predicting transport properties as 
envisaged here will be assessed by a comparison with numerical solutions of the underlying 
master equation for short chains containing four sites. The numerical effort needed for 
such solutions increases exponentially with system size (here expressed by the number of 
sites between the left and right reservoirs). For this reason, establishing the validity of an 
approximate solution is important for future applications. 

After defining our model in Sec. Ull we present in Sec. Illll a minimal description in terms 
of a 4-site model. Section |IV] introduces the TDFT method, while Sec. |V] examines its 
accuracy in comparison to "exact" numerical solutions. Also in Sec. |V] we present results 
concerning the performance of the system studied, focusing on the efficiency of its operation 
either as a motor working against a mechanical or electrical load, or as a chemical pump 
acting against a chemical bias. Section I VI I concludes. 

II. MODEL 

Our system is a 1-D "Fermionic" lattice gas with sites / = 1,...,M, time- dependent 
site energies ei{t), nearest neighbor hopping, and a nearest neighbor interaction V. By 
considering a "Fermionic" lattice gas, equivalent to local hard core repulsions, effects of 
saturation of site occupations are automatically included. The sites / = 1 and M can 
exchange particles with left and right reservoirs, L and R, that are sometimes represented 
below by indices and M + 1, respectively. We assume that these reservoirs exchange 
particles with the system with characteristic specified rates, but we disregard interactions 
between particles in the systems and those in the reservoirs, i.e. Vb,i = Vm,m+i = 0. By 
definition, reservoir particles are always in equilibrium and have fixed mean occupations. 



Pj 



(e-'^'^'^ + l) \ J = L,R (1) 



where fiL, ^r are the respective reservoir chemical potentials, and (3 = l/ksT. Rate equa- 
tions for the averaged site occupations p;(t); / = 1, ..., M; are given by 

-^ = {ji~i,i)t - (jM+i)t ' (2) 



where {ji^i+i)^ for / = 0, ..., M denotes the net average current from site / to /+1, to be derived 
from the underlying master equation. Note that (jo,i)t = {JL^)t ^^'^ {JM,M+i)t = {JM,R)t 
are currents from the left and to the right reservoir. In order to proceed, we need to specify 
the rates for configurational transitions consistent with the detailed balance condition. We 
adopt here symmetric rates 

w,j oc exp [f3{Ei - Ef)/2] , (3) 

where Ei and Ef is the total energy in the initial and final state, respectively. In par- 
ticular, 'bare' transition rates for elementary hops from / to / ± 1, that govern the single 
particle dynamics in the dilute limit, are given by ki^i±i(t) = iyii±iexp [f3{ei(t) — ei±i{t))/2]. 
Here i'i,i±i are frequency factors, which for simplicity are assumed independent of / for 
/ = 1, ■ ■ -M — 1 and represent the bulk frequency Ub, while at the system boundary we 
distinguish z/^ = z/o,i = ui^ and ur = vm,m+'^ = ^m+i,m from vb- Setting {vbY^ = 1 and 
/3 = 1 defines our units of time and energy. 

Certain special cases of this model deserve special attention. For static site energies 
El it describes aspects of passive transport, for example, through membrane channels.—"— 
Very recently, the nonlinear dc response and rectification in a single particle hopping system 
coupled to reservoirs were examined, including disorder effects.— Assuming pl 7^ Pr but 
taking ei independent of both time and space, one recovers a generalized asymmetric simple 
exclusion process (ASEP) model^S that contains a nearest neighbor coupling V (for reviews 
of the hard-core repulsion ASEP as well as the totally asymmetric simple exclusion process 



(TASEP) model, see Refs. |29| and|30| and references therein). 

Several driving modes, both stochastic and deterministic, that lead to unidirectional 
transport, were proposed in the literature.—*^ In this work we assume for simplicity a 
sinusoidal time dependence with frequency u = 2 tt/t, which should allow us to study 
transport efficiencies as a function of the typical time scale set by the modulation period 
r that characterizes the driving forces. With regard to the spatial asymmetry, common 
models are:— 

(a) Peristaltic or travelling wave-like behavior, where potential minima move in one di- 
rection, thereby dragging particles with them. 

(b) Sawtooth-like potential with oscillating amplitude ("flashing ratchet"). In such sys- 
tems particles are driven in the direction against the steeper potential slope. 



(c) Constant (/-independent) ac-force ("rocking ratchet") superimposing a non-symmetric 
static potential. Unidirectional transport directly results from steady state rectifica- 
tion properties of the static potential, as can be seen by considering the adiabatic limit 

Unidirectional transport in the latter model results from its DC rectification properties 
determined by the static potential, as can be seen by considering the adiabatic limit a; — )■ 0. 
The present work focuses on the first two mechanisms, (a) and (b). The model introduced 
above applies to multiply occupied channels driven by time dependent mechanical forces, 
and working as motors or as chemical pumps. As specific examples we consider the following 
situations: 

(i) Motor action: A constant load F is applied that changes the site energies relative to 
their intrinsic values {ei)F=o, 

Fl 
£; = (£/)f=o + p^— yy; / = 0,...,M + 1. (4) 

By this, the left and right reservoirs acquire a potential energy difference Er — El = F, 
while their chemical potentials are taken equal, p^ = p^. 

(ii) Chemical pump: In a pure chemical pump the chemical potentials in the two reser- 
voirs are different, e.g., fiR > fii {pr > Pl)-, while the mechanical load F vanishes. 
Mathematically this case differs from the model dl]) by the fact that, unlike the linear 
potential change in (|1]), the chemical potential difference /i/j — ^l will in general not 
give rise to a linear distribution of local chemical potential changes along the channel. 

In both cases we can discuss different measures of machine performance and optimization 
schemes. First the output of the machine operation can be obtained from the average current 

1 r 1 ^^ 

^av=-/ dtj^^Y.^3lHl)t (5) 

where, at steady state, the average (over realizations and over a modulation period r) 
currents between neighboring positions do not depend on position. The useful work output 
is then 

Wont = J..F. (6) 



for the pure motor action, and 

Wont = Ja.vifJ'R " /"l) • (7) 

for the pure chemical pump. Note that combined effects of mechanical and chemical biases 
can be considered, in which case H^out = Ja.vifJ'R — A^i) is the sum of ([6]) and (171), the 
electrochemical potentials /i| j^ being defined according to (140|) . In all cases the work input 
can be calculated from 

M 



dei{t) 

'"'it 

This yields the conventional "efficiency"— 



'^-=7r*E^«w («) 



Wont .^. 

ri = — . (9) 

W, ^ ^ 



m 



Aiming for maximum efficiency is one criterion for optimization. Alternatively, irrespec- 
tive of the amount of input energy, one can ask for the maximum F, where Jav changes 
sign ("reversal potential"), or for the maximum current Jav > against the load F > 0, 
corresponding to the maximum rate of transfer of particles from L to R. 

III. FOUR-SITE MODEL 

Setting M = 2, we have a 4-site system consisting of a channel with sites / = 1, 2 in 
contact with boundary sites L and R. In this minimal ratchet model open to reservoirs the 
system is driven by modulating the site energies ei{t) and e2{t). An incipient peristaltic 
modulation of site energies can be realized by a phase lag in the oscillation of B2{t) relative 
to ei{t), 

eL = (10) 

e, (t) = ef ) +A[1 + sin(wt)] + F/3 (11) 

e2(t) = ef +A[1 + sm{cot - n/2)] + 2F/3 (12) 

eR = F. (13) 

Here, the energies e\ represent a constant energy shift of the channel's interior relative to 
the reservoirs. On the other hand, the in-phase oscillation where Eqs. (ITT|) and (IT2l) are 
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replaced by 

£^{t) = ef^ + 2A[1 + sm{ujt)] + F/3 (14) 

£^{t) = 4°^ +A[1 + sm{ujt)] + 2F/3 (15) 

(keeping fITU]) and (IT5]) ). corresponds to a "flashing" ratchet of the type of a discrete "saw- 
tooth" potential. Both of these driving modes favor a current to the right. 

Unidirectional flow induced by such driving schemes can be most simply investigated in 
the independent particle model where the average site occupations Pi{t) and P2{t) evolve 
according to the linear rate equations 

-TT = kL,i{t)pL - ki^L{t)pi + k2,i{t)p2 - ki^2{t)pi (16) 

-7^ = kR^2{t)PR - k2,R{t)p2 - k2,l{t)p2 + ^1,2(^)^1 (17) 

with prescribed populations pl,r at the boundary sites. A more realistic model should take 
into account interparticle interactions. Here we consider both hard core (site exclusion) and 
nearest neighbor interaction, denoted V. 

The following section (Sec. IIVI) describes an approximate approach to the kinetics of 
such models based on the (classical) time-dependent density functional theory. An exact 
approach, feasible for the present small system is based on rates equations written in the 
system states representation. A system state (721,722) is defined in terms of the occupations 
of sites / = 1 and / = 2. Obviously our system is fully characterized by the four states (0, 0), 
(1,0), (0, 1) and (1,1), with the corresponding energies 

^io(t)=£i(t) (18) 

Eoi{t)=e2{t) (19) 

En{t)=ei{t)+e2{t) + V. (20) 

Eqq may conveniently be set to 0. We assume that transitions between these states proceed 
only by single particle steps, so that no direct transition takes place between states (0, 0) 
and (1,1). It should be emphasized that in using these energies to determine rates one 
needs to take into account the change of energy in the reservoir. Thus, the total energy 
change AE for the transition (0,0) — t- (1,0) is Eiq — Sl = Ei—El and AE for the transition 
(1, 0) — 7- (1, 1) is Ell — EiQ — Er = E2 + V — Er. The kinetic equations for the average 



population of these system states are 

dP (f) 

-^^ = KooMt)Poo{t) + Koi,io{t)Poi{t) + i^ii,io(t)Pii(t) (21) 

- (Kio.ooW + Kiofiiit) + ^^10,11 W) Pioit) 

^^^ = Koo,oi{t)Poo{t) + K,o,oiit)Pioit) + /^ii,oi(t)Pii(t) (22) 

- {Koi,oo{t) + Koi,io{t) + KoiAiit)) Poi{t) 

dP (t) 

-^ = Kio,ii(t)Pio(t) + i^oi,ii(t)Poi(t) - (/^ii,io(t) + /^ii,oi(t)) Pii(t) (23) 

and normalization implies that 

Poo = l-(Pio + Poi + Pii). (24) 

Using (124|) in ( l2Ti) -( l23l) yields three inhomogeneous equations with rates that are readily 
obtained from ^ and ([IS]) -([20]). Clearly, Km^wif) = k2,i{t) and Kio,oi(i) = ki,2{t). The 
remaining ones involve the bath densities, for example, 

KioMt) = i^Lexp[/3(ei(t) - Sl) /2]{1 - pl) (25) 

Koo,io{t) = vl exp[-/3(£i(t) - ElM'Avl (26) 



Once the solutions to Eqs. ([2Tl) - ([2ll) have been found, the average populations of individual 
sites are obtained from: 

Pi(t) = Pio(t) + Pn(t) (27) 

P2(t) = Poi(t) + Pii(t). (28) 

and the currents Jl(^) = (jl,!)^; Jnify = (J2,/?)t can be calculated from, e.g. 

Jiit) = -f^oo,ioPoo ~ -^io,ooPio + -^^'oi.iiPoi ~ -^ii,oiPii) (29) 

Jnit) = A'oi,ooPoi - AVoiPoo + i^ii.ioPii - AViiPio- (30) 

Note that results based on Eqs. ( [2T]) -( [24]) differ from those of the independent particle 
model, Eqs. ( [T6[) - ([T7j) . even in the limit V = 0, because unlike the latter they incorporate 
hard core interactions. The behavior based on Eqs. ([T6]) - ([T7)) is expected only in the highly 
dilute limit when particle encounters on the same site are negligible. 
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IV. TIME-DEPENDENT DENSITY FUNCTIONAL THEORY (TDFT) 

TDFT is a local equilibrium approximation, in which the non-equilibrium character of 
the distribution function is manifested in the space- and time-dependent local fields acting 
on the single particle density. Given a distribution function of this type, density functional 
theory assumes that correlators determining the currents (j;,/+i)^ in Eq. (|2]) are functionals 
of the single particle density, and moreover, that this functional dependence is the same 
as in the equilibrium case. This allows us to express the mean currents in ([2]) in terms of 
the densities Pi{t) and therefore to arrive at a closed system of nonlinear rate equations. In 
the present problem this latter step can be carried through exactly because the free energy 
functional for 1-dimensional lattice gases with short range interactions is known.— i^ 

The rate equations ([2]) for the averaged site occupations pi{t)] / = 1, 2; in the 4-site 
model reduce to 



dp 

It 

dp2 

dt 



-^ = {JL,l)t - {jl,2)t 
= {jl,2)t - {j2,R)t 



(31) 
(32) 



From the general formulation given in Appendix |A] for channels with arbitrary length M 
in the presence of symmetric rates ([3]), we obtain by specializing to the present model with 
M = 2, 



OL,l),= (l-pi+i^pg 



~kL,ie^>''-~h,L 



(2)' 
(4) , 



(3) 

k2J I %y 1 - ~kR,2 e^''^ 



(33) 
(34) 
(35) 



with K = VC- 1; C 



-I3V 



(for details see appendix |A]). For the sake of simplified notation. 



time arguments on the right hand side of these and subsequent equations are suppressed. 
The nearest neighbor correlators pg i i^ these equations are defined by Eq. flA16p with / = 1. 
Note that they are directly related to the average occupations of system states. Indeed, 
Pii = P21, Poi = P2 i5 -Pio = P2 1 ^^d Poo = P2 1- Following the Appendix, the two- point 
correlator P2 1 — ('^2'^i) is explicitly found as 



,(1 



mt) 



2(1-0 



iPi+P2) (1 - c) - 1 + \/[ipi+P2) (1 - c) - 1]' + 4piP2c (1 - 



(36) 

Clearly, setting all currents (155]) to fl5H]) equal to zero implies the equilibrium condition 
£l + f^L = £r + fJ'R between both reservoirs. This is seen most directly from the underlying 
equation (lASp . 

In the limit V — > one recovers P21 -^ PiP2- Then Eqs. ( 133|) -( 135|) become 

(jL,i)t = kL,iit)pL (I-Piit)) - h^Lit)piit) (I-Pl) (37) 

{jia)t = kl,2{t)Pl{t) {l-p2{t))-h,l{t)p2{t) {1-Pl{t)) (38) 

(j2,R)t = k2,R{t)P2{t) (I-Pr) - kR,2{t)pR (1 -P2(t)) (39) 

showing that site blocking effects are incorporated in a mean-field like manner. For general 
V, effective blocking factors for sites / = 1,2 can be read directly from (l33l) and (!35ll . 
Considering site 1, the effective blocking factor that appears in the current (JL,i)t is l—pi + 
K P2 j, which reduces to 1 — pi for l^ = and to 1 — (pi +P2) for ^ — ;■ 00, because P2 1 "^ 
in the latter case. Hence, for V — !■ 00 the occupation of the nearest neighbor site, p2, enters 
into the effective blocking factor additively with pi. The last conclusion follows from ( l36ll 
provided that pi+ P2 < 1, which obviously holds for V -> 00 provided that fiR and fiL stay 
finite.— By this, and based on the identification Pn = P21, it is straightforward to show 
that for y — > 00 the expressions for the currents (l33l) - (!35|) and hence all results from TDFT 
become equivalent with those of Sec. IIIII in the same limit (i.e. for Pn — )■ 0). A similar 
conclusion holds for V — ?■ —00, where -R'lo,!! — ;■ C)0. According to ( I23l) this limit requires 
that Pio —7- 0. Similarly, Pqi — ^ 0, and Pn = pi = p2 = 1 independent of time is seen to 
be a consistent solution. In this limit all currents are zero as the channel is jammed by the 
presence of two particles. On the other hand, (136|) predicts that P21 ^ max(pi,p2), which 
becomes unity since pi,p2 — ^ 1- It follows that the TDFT correctly yields zero currents in 
this limit. 

The following section provides a quantitative test of TDFT against the exact results ob- 
tained from Eqs. (12T]) -( |2^ . When applied to the 4-site model, this approximate method is 
not simpler than the exact treatment, however, its advantage lies in the immediate applica- 
bility to chains of arbitrary length M.— 
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V. RESULTS 

A. Peristaltic driving, dilute limit 

Now we investigate a 1-D open channel under peristaltic driving by numerically solving 
Eqs. (inD-dll in Sec. [ml and Eqs. (13ID-(I36]) in Sec. IV] with the input Eqs. (Il0])-([l3]). We 
first consider the dilute limit, approached by setting /ii = /i_R = —5, i.e., pl = Pn. ~ 0.0067. 
Calculated efficiencies as a function of the mechanical load F show a maximum before they 
drop to zero and to negative values (current reversal), see Fig. [TJ Clearly, for strong dilution, 
site blocking effects and the interaction V become irrelevant. Hence, the efficiency curves 
in Fig. [1] become indistinguishable from the predictions of the independent-particle model 
Eqs. ( !T6|) and (TT7I) . Also in this dilute limit, the TDFT results in Sec. lIVI perfectly agree with 
the numerical solutions of the exact equations fl2T|) - fl2^ . Indeed, in this limit Pn becomes 
negligibly small, and, as one can easily verify, both the exact equations (1^T]) - (I2^ and the 
TDFT equations f l^ - fl5^ become fully equivalent to the linear rate equations flT^ and 



Figure [T^ also reveals a significant enhancement of the efficiency for a channel which is 
made more attractive relative to the reservoirs, by changing the static part of site energies 
in (ITTl) from e]^ 2 = to e^ 2 = —2.0. The value ^i 2 = —2.0 is therefore used in all our 
subsequent calculations with peristaltic driving. Apparently, attracting more particles to 
the channel interior overcompensates for the effect of a less favorable exit rate k2fi. 

Direct insight into the peristaltic mechanism is gained from Fig. [T|d. During a time 
window where Si{t) ^ 1 particles from the left reservoir have access to the channel so that 
the left current </l(^) and pi{t) increase with time. (This happens for 1.5 ^ t ^ 2 in the 
plot of Fig. [!}).) During the subsequent upward movement of level 1 density from this level 
flows both back to the reservoir, rendering Jiit) < 0, and in the forward direction to level 2 
, which is lower in energy because of the phase lag 7r/2 between both levels. Therefore P2{t) 
increases, but as long as e2{t) ^ 1 there is also a flow to level 2 from the right reservoir, 
reflected by the negative peak in the right current JR{t). This influx from the right is the 
reason why the p2-peak actually gets higher than the pi-peak. Upward movement of level 
2 in turn causes the subsequent positive peak in Jnit), while a decay of p2(^) to the left is 
prohibited as long as e2{t) < ei{t). The peaks in the time-dependent densities and currents 
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can be shown to be more pronounced for e^^ 2 = ~2 in comparison with £^2 = 0; leading to 
the increased efficiency displayed in Fig. [1^ for the more attractive channel. Since during 
the upward movement of level 1 its population decays both towards the left and right, we 
expect that this mechanism can translocate a particle with probability not larger than about 
0.5 within one cycle. 



B. Peristaltic driving, moderately dense system, pl = Pr 

Figure Wp shows that at low densities, a large amplitude peristaltic oscillation Eqs. ( TTOj) - 
( IT3I) tends to locate particles on that site which momentarily has lower energy. This is evident 
from the figure as the overlap of the two peaks for pi(t) and P2(t) is small. Figs. [2] and |3] show 
results obtained at higher densities, imposed by setting p^ = p^ = 0.5 (/i^ = /ir = 0), where 
effects of correlations between transporting particles are expected. At these densities the 
overlap between the pi{t) and P2it) peaks is larger, as seen in the case l^ = in Figs. Et and 
|3b , which show slightly broader peaks and overlap regions than in Fig. [T]d. Consequently, 
many-body effects become important, as discussed below. 

It may be intuitively expected that simultaneous presence of particles on both sites 1 
and 2 should cause reduction of pump efficiency, because these particles block each other, 
impairing the peristaltic driving. Hence one expects that efficiency will increase as V changes 
from attractive (V < 0) to repulsive (V^ > 0), where double occupation of the channel is 
increasingly suppressed. These arguments are supported by calculations of rj and Jav for 
non-dilute systems. Figure [2] shows results for Pl^r = 0.5, obtained from both the exact 
method described in Sec. IIIII and the TDFT approximation. Sec. IIVI Before addressing the 
quality of the TDFT approximation, let us focus on the main features in Fig. [2] common to 
both treatments. 

Comparing the V = case^ shown in the inset to Fig. |2^ with the independent particle 
results (dilute limit. Fig. [1^) we indeed observe that mere site blocking leads to a decrease 
in rj. This trend in rj is enhanced for increasingly negative V (attractive interparticle inter- 
actions) but is reversed with increasingly positive (repulsive) V, again in agreement with 
the above expectation. Interestingly, for \^ = 10, 77 significantly exceeds the result for the 
dilute limit discussed before, see Fig. [1^, by about 30 percent. One should note that a 
strong repulsion as in the case V^ = 10 is practically equivalent to the condition that the 
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two-site channel can only be singly occupied or vacant. This "single particle limit" appears 
to optimize active transport in our model. 

Similar conclusions hold with respect to the average current Jav, plotted against F in 
Fig. [2]3. The current increases for increasingly repulsive (positive) V and becomes negligibly 
small for increasingly attractive interaction. For ^ = 10 the reversal potential (the load F 
for which Jav vanishes) is larger by about 15 percent than the corresponding value in the 
case y = 0. For strong attractive interaction, e.g. V = —5 , see Fig. [2], both current and 
efficiency are strongly damped. 

Regarding the comparison of both methods, it is evident from Fig. [2] that the TDFT can 
well account for important qualitative features in our model of a driven open channel. At 
V = 0, the TDFT becomes identical to ordinary mean-field theory because all correlators 
factorize, leading to Eqs. ( J37l) - (j39l) . One should note that quite generally non-trivial and even 
long range correlations can be induced in systems when driven far away from equilibriunt^'^. 
Such correlations are ignored in the TDFT approach. From the inset to Fig. [2], we see that 
for V^ = TDFT overestimates the maximum efficiency r] by about 14 percent relative to 
the exact calculation. This error shrinks when V increases from zero, because the TDFT 
accounts for the concomitant suppression of double occupation of the channel. For strong 
coupling, V ':^ 1, the TDFT applied to our model asymptotically becomes exact because, 
as discussed before in Sec. IIVI it correctly describes the one-particle limit. 

For practical applications, attractive interactions are less interesting than repulsive ones 
since they tend to immobilize a pair of particles on sites 1 and 2, thus diminishing the 
current. For the sake of comparison of the two methods, however, we have included in Fig. [2] 
examples with V^ < 0. An almost quantitative accuracy of the TDFT is observed when 
the attraction becomes as strong as V = —5, again in accord with the previous analytical 
arguments in Sec. IIVI 

The development of correlations inside the channel of the type discussed above is illus- 
trated in Figs. [St; and[3li, which show the difference -Pii(t) — Piit)p2{t) within one driving 
period. The load is taken as F = 1.5, where the corresponding efficiencies are close to their 
maximum. For V^ = 0, see Fig. [3|i, peaks of either sign appear in the exact results in parallel 
with the peaks in pi{t) and P2{t), but they are narrower than the latter. On average, these 
current- induced nonequilibrium correlations are positive, i.e., attractive. However, already 
under the mild repulsive interaction, V = 1, the negative peak takes over. For strong re- 
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pulsion, Pii practically vanishes, so that the curve for V^ = 10 simply reflects the negative 
product of the densities pi and p2- Obviously, the TDFT reproduces such correlations better 
for larger V. 

The influence of the driving period r (in fact t/ub as discussed above) is displayed in 
Fig m Unless otherwise specified, here and in the following we only show the exact results 
from Sec. Illli In all cases studied, however, we have verified that the TDFT performs with 
similar accuracy as for the foregoing plots. Obviously, except for the driving mode itself, the 
driving period is a fundamental parameter that determines the pump performance. For our 
model both the the efficiency t] and the current Jav go through a maximum as a function of 
T. To maximize rj, a smaller load requires slower driving, as seen in Fig. H^. Increasing F, 
the maxima in rj shift to shorter r. On the other hand, from the point of view of maximizing 
Jav, even shorter r are required, e. g. r ~ 1.2 for F = 2 (see Fig. |Dd). This illustrates that 
optimization of pump performance in general depends on the optimization criterion for a 
particular application. The choice r = 2 in most of our calculations is found to be a good 
compromise between these different optimization criteria. Yet another criterion might be to 
maximize the number of particles transmitted within one cycle, which is r Jav Figure Sb 
suggests that in our model this number will be limited by about 0.5, as argued before. It 
should also be noted that r-dependent efficiencies and the peak structure in the Jav versus 



r curves qualitatively agree with Refs. 
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and 



16l . where a model based on repulsive on-site 



interactions and periodic boundary conditions was considered. 

C. Peristaltic pump; effect of chemical driving, pl 7^ Pr 

It is well-known that passive transport in the linear response regime at constant temper- 
ature is governed solely by the gradient in the external electrochemical potential 

fi' = F + fi. (40) 

Far from equilibrium, however, the mechanical and chemical potential F and /i play separate 
roles in determining the kinetics. This difference stems from the different ways at which 
these imposed biases are expressed inside the channel. The "mechanical bias" F = er — Sl 
is assumed to fall linearly along the channel, as described by Eqs. ( TT0|) - (TT5|) . while the 
chemical bias A/x = (j^r — fJ.L is assumed to affect only the external reservoirs R and L. 

14 



Therefore, the average current Jav in our model depends differently on these two external 
variables, Jav = Jav(-^, A/i). This is illustrated in Fig. |5l where the current is plotted 
versus the difference A/i^ = fi'j^ — fi^ in the electrochemical potentials of the two reservoirs 
under different partitionings of A/i^ with respect to the mechanical and chemical load. In 
these calculations pl = 0.1 is fixed {fii ~ —2.197), whereas pr > pl is variable upon 
varying A/i. Shown are only the two limiting cases Jav(0,A/i) and Jav(-F, 0), labelled as 
F = and Afi = 0, respectively. Note that in biophysical systems pumping of charged ions 
across a membrane will generate a voltage drop so that in principle the full two-variable 
characteristics Ja,^{F, Afi) will enter. As seen from Fig. Eti', the current Jav(0, A/i) for small 
bias is less sensitive to changes in Afi in comparison with the F-dependence of Ji^^{F,0). 
This is because, after time averaging, pft enters the net current from site 2 to i? mainly 
through the blocking factor 1 — p/j in equation fl25|) . implying only a mild p/j- dependence. 
However, when pji > 0.5 the current Jav(0, A/i) drops more steeply towards negative values 
than Jav(-^, 0). In contrast, as in Fig. [HI the drop in Jav(-^, 0) down to the reversal potential 
is not far from linear. The corresponding TDFT results were found to be in excellent 
agreement with these results. Note that the small concentration pl = 0.1 implies that the 
interaction V has only a minor influence on Jav(-^; 0). 

The difference in the channel's response under a chemical versus mechanical load also 
becomes apparent when we compare the respective efficiencies, which are plotted in Fig. [6] 
against A/i"^ = At^ — /i|. As for Fig. [5] we have taken p^ = 0.1 fixed and pn variable. The 
chemical efficiency, defined by Eqs. ([9]) and (0), is displayed by the curves labelled F = 0. 
Their maxima and the corresponding chemical reversal potentials are systematically lower 
than those referring to the mechanical efficiency (A/i = 0, F ^ 0). 

D. Flashing ratchet potential 

An important question is to what extent the previous results for peristaltic driving will 
change when the drive mode changes. We exemplify this for the flashing ratchet potential 
Eq. ( IT^ . Here we use the parameters A = 2.0 and £12 = —1-0 that correspond to a maximum 
in e2{t) of the same height relative to the bath levels as for peristaltic driving with A = 5.0 
and ^1 2 = —2.0. Figure [7] shows F-dependent mechanical efficiencies 77 and currents Jav 
The qualitative appearance of the curves including their dependence on V is analogous to 
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Fig. [21 but the absolute performance is considerably less good than for peristaltic driving. 
The efficiency rj and the average current Jav are seen to be lower by about an order of 
magnitude and a factor of 5, respectively, in the flashing driving mode. Thereby we have 
verifled that choosing r = 2, the efficiency curves ?7(-F) in Fig. [7] are near optimum. We 
also computed the analog to Fig. \5^ for the flashing ratchet, see Fig. \5jp. Currents again 
are much smaller than in Fig. [5^. A notable feature is the near agreement of the curves 
for chemical {F = 0) and mechanical (A/i = 0) bias, presumably because of a fairly regular 
distribution of induced local chemical potential changes along the channel in the case of the 
sawtooth potential. 

VI. SUMMARY AND CONCLUSIONS 

Active transport was studied within a stochastic model for an open two-site channel, 
driven by time- dependent site energies and coupled to left and right reservoirs. Exact rate 
equations for this problem were established and solved numerically. Despite the simplicity 
of our model, it allowed us to examine important conditions for efficient transport, which 
we expect to be of relevance regarding both the design of synthetic active channel devices 
and, potentially, the functioning of biological motors and ion pumps. 

The main findings which emerge from our investigations are summarized as follows: 

i) Because of the coupling to reservoirs with prescribed chemical potentials, the number 
of particles inside the channel fluctuates, allowing multiple channel occupation. As 
a consequence, transport depends on interparticle interactions, described by a local 
hard core repulsion and a nearest neighbor interaction constant V. l^-dependent motor 
and pump efficiencies are studied systematically, showing that they become optimized 
when multiple occupation of the channel gets suppressed by a large repulsive V. As 
intuitively expected, this effect saturates and the efficiencies become independent of 
interparticle interaction when l^ — )■ oo. In particular, the V^ = 10 results displayed 
in Fig. [2] were found to represent this saturation limit. On the other hand, choosing 
an even higher modulation amplitude {A = 10 instead oi A = 5 used in the foregoing 
calculations for peristaltic driving), the maximum efficiency was found to increase only 
slightly, whereas the reversal potential increases considerably. 
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ii) Different performance criteria can be formulated, based, for example, on an economical 
use of the input energy in producing mechanical or chemical work (efficiencies 77), or 
on maximizing the output "power", i. e. the current. These measures of efficiency, as 
defined in this paper are not intrinsic properties of the system, but obviously depend 
on the imposed load. It is also possible to define an intrinsic measure, 77^, of an "ideal" 
efficiency by (for a mechanical bias; a similar definition can be used in the chemical 
bias case) 

_ J..(F^O)F(J^O) 

where Jsv{F = 0) is the flux at zero bias while F{J = 0) is the mechanical bias 
for which the current vanishes. The product Jav(-^ = 0) F{J = 0) is the analog of 
the product of the short-circuit current and the open-circuit voltage that provides an 
upper bound to the useful work that can be extracted from a given voltage source. 
The choice of Win_{F = 0) as the denominator in Eq. (14T|) is to some extent arbitrary, 
since the input work in our model is defined through the energy absorbed by the 
system (in contrast to situations encountered, e.g., in photovoltaic devices, where the 
input work is defined by the incident rather than the absorbed energy). Other choices, 
e.g. using Wi^iJ = 0) as denominator in ( HTl) could be made. As shown in Fig. [8] 
these different definitions yield different results for what may be regarded as the ideal 
machine performance. In either case, the performance is sensitive to the interparticle 
coupling and to the driving period r, and different efficiency criteria require different 
r for optimization. 

iii) In the far-from-equilibrium processes considered here, generated currents depend on 
the partitioning of the electrochemical potential difference A/x*^'^^ = F + A/i between 
the reservoirs in terms of the mechanical and chemical load F and A/x. In other words, 
these biasing attributes play different roles in the kinetics. 

iv) Peristaltic driving reflected by a phase shift in the oscillating site energies yields much 
better performances than an oscillating sawtooth- like potential (flashing ratchet). 

Some of these findings, in particular the possibility of efficiency enhancement by repulsive 
interactions and the above-mentioned dependencies on the driving period agree on a qual- 
itative level with earlier work,— ^^^ despite considerable differences in the model design and 
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in details of the results. We take this as an indication that the corresponding conclusions 
formulated above should hold rather generally. 

It would certainly be interesting to extend these investigations to longer channels, with 
the aim to analyze collective effects in overdamped systems under time- dependent driving. 
The full problem based on the master equation, however, soon becomes intractable as the 
length M increases. Results based on the time-dependent density functional theory, TDFT, 
nicely agree with the exact solutions for M = 2, suggesting that this approximation could 
yield reliable results also for arbitrary channel lengths. The underlying equations of motion 
for arbitrary M are presented in the Appendix and will be evaluated in forthcoming work. 

It should be emphasized that, unlike equilibrium properties, the behavior of the systems in 
non-equilibrium steady states depends, sometimes sensitively, on the assumed kinetic model. 
We have repeated some of the calculations for two more kinetic models, a normalized kinetic 
model, where the rates between two levels with energy separation AE > are taken to be 

(42) 
yuowny ea"'^ -|- e~^^^ 

and an asymmetric model where 

6 - -.- 

fc(up) = ^— ; k{down) = ^— . (43) 

1 + e-2^-^ 1 + e"2^^ 

Because rates based on Eqs. fj42|) or fHSl) are smaller than those from Eq. ([3]) longer mod- 
ulation periods are needed in order to get similar behavior. Apart from this rescaling of r 
the qualitative behavior of these models remains as with our original model, Eq. ([3]). 

Some other important questions are deferred to future work. In addition to the average 
currents studied in the present paper, current fluctuations are observed in experimental 
studies of channel transport, and their theoretical evaluation is of interest. In particular, 
generalizing available methodologies for carrier counting statistics to models with time- 
dependent driving provides an interesting direction for future investigation. 
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Appendix A: Currents from the TDFT approach 

In this Appendix we recall the main steps of the TDFT for "Fermionic" lattice gases^^i"— . 
At the same time we generalize the derivation to systems with arbitrary time dependent site 
energies e^ and kinetic properties governed by a class of transition rates, where 

Wi^k (n) = Vi^kTiiVkg {Hi, Ht) (Al) 

represents the rate for a particle hop from site i to a nearest neighbor site k. Occupation 
numbers are denoted by rti with ri; = or 1, vi = 1 — ni and n = {rii}. Vi^k are frequency 
factors. The function g depends on the energies Hi and H^ of the particle before and after 
the hop, respectively. Focusing on this transition between sites i and k, the total lattice gas 
Hamiltonian is decomposed as 

H (n) = Hi^k + riiHi + UkHk + riiUkVi^k , (A2) 

where Vi^^s are pair interaction parameters. Hi^k denotes the total energy omitting contri- 
butions from sites i and k, while 

Hi = Si- fitot + ^ Vi^ini . (A3) 

Here fitot is some reference chemical potential. 

Clearly, Hik, Hi and Hk depend only on occupational configurations n excluding sites i 
and k. For the following it is convenient to introduce the quantity 

G,,k (n) = g {H,, Hk) exp {-(3H,) . (A4) 
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which by the detailed balance condition is required to be symmetric in its indices, 

G,,fc (n) = Gk, (n) • (A5) 

If transition rates depend only on the particle's initial energy,— then Gik (fi) = 1. On 
the other hand, calculations in this work are based on the symmetric choice g {Hi, H^) = 
exp[l3{Hi-Hk)/2]oT 

G,,fc (n) = exp [-/3 {H, + H^) /2] (A6) 

The TDFT is based on a local equilibrium approximation for the distribution function. 



Pi- (n, t) = -^ exp [-^H^^ (n, t)] ; H^^ (n, t) = if (n) + V h, (t) m . (A7) 

Its non-equilibrium character is represented in terms of time-dependent single-particle fields 
hi (t) that are associated with the time dependent local chemical potential, see flAip . In 
view of Eqs. (IAip - (lA5p and using nf = nf, UiVi = 0, the average current from i to k can be 
expressed as 

{ji,k)t = ^i,k J^ [niVkQ {Hi, Hk) - ViUkQ {Hk, H;)] P^"" (n, t) 

n 

A,k J^ Gi^k (n) [riiVk - ViUk] exp [-(3 {hiUi + hkUk)] exp {-(3Hf^ (n)) 



n 



where H^^ is to the first term in ( JA2I) when applied to the corresponding decomposition of 
H'^^. Summation over rii and Uk yields 



{j,^k)t = M,^k {t) [A, (t) - Ak (t)] (A8) 



where 



and 



A, (t) = exp [-(3 {h, {t) - /ieq)] (A9) 



^^-^ ^^) = 1^ E G.,k (n) e-^^^l^'^^e-^^^ 
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The summation over fi can be continued to a summation over all n after including vacancy 
occupation numbers for sites i and k . Thus 

M,,fc (t) = i/,,fc {viVkG^^k (n)), e-^>'"" , (AlO) 



which is symmetric, Mi^k{t) = Mk,i{t), in view of (IA4p . Note that ^tot cancels in (lAlOp . 



Density functional theory is employed when computing average densities pi (t) = (nj)^ from 
(IA7p . The result is the "structure equation" 

ei + hi + ij,i{p) = i^tot (All) 

which allows us to eliminate hi in favor of the local chemical potential /ij (p) = |^; F (p) 
being the free energy functional associated with H (n), with p = {pi} (for details see Ref.llSl). 
It follows that 

A, (t) = exp [/3 (Ei + ^i^)] , (A12) 

which can be interpreted as local activity. The essence of Eq. flA8l) now becomes obvious, 
namely a factorization of the average current into a thermodynamic factor (difference of 
local activities) and the kinetic coefficient flAlOp . which is symmetric, Mj ^ = M^ j. A closed 
system of equations determining pi (t) is obtained by combining Eqs. (lA8p - (lA12p with the 
equation of continuity 

^ + Efe.>, = 0. (A13) 

k 

Now we apply that method to the model of Sec. ([111), which is a finite one-dimensional 
channel of M sites with nearest neighbor interactions, coupled to reservoirs L and R. From 
(IA6D and (lATOJ) we obtain 



Mi^i+i (t) = ui^i+i exp [-13 (ei + ei+i) /2] {vivi+i exp [-/3 (Vi-i^mi-i + ¥1+1^1+2^+2)]) ^ (A14) 

where Vi-i^i = V for / = 2, ...M. Note that the reservoirs were taken to have no interaction 
with the channel, Vb,i = Vl,i = and Vm,m+i = Vm,r = 0. Moreover, it exactly holds thai^^ 

Pfii (p) = In + In + In (A15) 

^ P^ PiPi+i,i PiPiJ-i 
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(n) 



The correlators pl^\ ^ are defined by 



p!+i,/ = {vi+iui) =pi- pS:J.\ ; (A16) 

p!+i,z = {ni+ivi) = pi+i - p^,^\i (A17) 

vf+i,i = (vi+ivi) = l-pi- pi+i + p{+\^; (A18) 

with ri; = or 1 the occupation number for site / and vi = 1 — ni is vacancy occupation 
number. These quantities pi^\ ^ can be determinedi^ from the "quasichemical condition" 
Pi+iiPi+ii ~ Pi+iiPi+ii e"''^ by solving a quadratic for, e. g., Pij^n = {rii+ini). Equation 
( IA15I) is written as a sum of three terms such that the first term corresponds to a non- 
interacting lattice gas, while the second (third) term is the contribution of the right (left) 
neighbor to the chemical potential at site i. 

Equations (1A14P - (]A15|) together with (]A12p determine the currents (JASP in the interior 



of the system. Evidently, flA14p involves higher order correlators, actually up to 4-point 
correlators, because of exp (an) = 1 + n (expa — 1) for n = or 1. Knowing -F(p), these 
correlators can be expressed^^ as functionals of p. 

The average currents from and to the reservoirs require special attention. First, because 
interactions between system and reservoirs are disregarded, the chemical potentials for the 
outermost sites / = 1 and I = M oi the channel satisfy 

(2) (3) 

/3/ii(p) =ln^; /3/iM (p) = In ^^p^ (A19) 

P2,l Pm,M-1 



This follows from the fact that all correlators in flAlSP involving reservoir sites factorize so 



that in ( 1A15I) contributions of the interactions between sites I = 1,M and the reservoirs 
vanish. Secondly, the kinetic coefficients that enter the currents between system and are 
obtained from (using Eq. (IA14p . 



Ml,i (t) = ul exp 



_^.e, + e^ 



(t;z.)(t'iexp(-/3^)) , (A20) 



Using again exp (an) = 1 + n (exp a — 1), the last factor in (1A20I) simplifies, 

v^ exp (-/3^) ) = 1-Pi + Kp'il (A21) 



with K = exp (— /3V"/2) — 1. An expression analogous to flA20p is obtained for MM,R{t). 
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After insertion of ( 1XT9|) -( 1MT]) into (IXT2D and ([MD we obtain 



iJL,l)t 



i-p^ + Kp'ii 






(2)' 
P2,l 



,^2 



(4) 



(A22) 



where {vl) has been absorbed in the attempt frequency entering ki^i and ki^L by setting 

^L = {vl) ^l- 
Similarly, 



(2) 



{3M,R)t =[1-PM + Kp\'_^ 



(3) 



(A23) 



with vn = (vfi) vr in the definition of A;_r,m and A;m,_r- Note that apart from modified attempt 
frequencies i'L,Ri the only reservoir properties entering the theory are the chemical potentials 
^L,R- In this way and by using flAlSp . we end up with a closed system of equations for an 
open channel of arbitrary length, coupled to reservoirs. 

Specializing to the 4-site model (M = 2) of Sec. Illlt the above expressions for the reservoir 
currents coincide with Eqs. ( I33l) and (l35l) in the main text. Moreover, to get the current 
{ii,2)i inside the two-site channel, we use ( 1A14I) to obtain 

'(^1 + ^2) 



Mi^2 if) = ^ exp 



-/3 



(^'l^2)i 



(A24) 



.(4) 



with {v\V2i^ = ^2 1- Together with (IA19p . this leads to ([3 



REFERENCES 

^F. Jiilicher, A. Ajdari, and J. Prost, Rev. Mod. Phys. 69, 1269 (1997). 

^P. Lauger, Electrogenic Ion Pumps (Sinauer Associates, Sunderland, MA, 1991). 

^B. Hille, Ionic channels of excitable membranes (Third edition, Sinauer Associates, Sun- 
derland, MA, 2001). 

^E. Muneyuki and T. A. Fukami, Biophys. J. 78, 1166 (2000). 

^V. Kaila, M. Verkhovsky, G. Hummer, and M. Wikstrom, Proc. Natl. Acad. Sci. 105, 6255 
(2008). 

^J. Rousselet, L. Salome, A. Ajdari and J. Prost, Nature 370, 446 (1994). 

^Z. Siwy and A. Fulihski, Phys. Rev. Lett. 89, 198103 (2002). 

^M. Galperin, M. A. Ratner and A. Nitzan, Nanoletters 5, 125 (2005). 



23 



^B. L. Altshuler and L. I. Glazman, Science 283, 1864 (1999). 
i°P. Reimann, Phys. Rep. 361, 57 (2002). 

"P. Haenggi and F. Marchsoni, Rev. Mod. Phys. 81, 387 (2009). 
^^I. Derenyi and T. Vicsek, Phys. Rev. Lett. 75, 374 (1995). 
^^S. Savel'ev, F. Marchesoni and F. Nori, Phys. Rev. E 71, 011107 (2005). 
^^P. Reimann, R. Kawai, C Van den Broeck and P. Haenggi, Europhys. Lett. 45, 545 (1999). 
i^F. Slanina, Europhys. Lett. 84, 50009 (2008). 
i^F. Slanina, J. Stat. Phys. 135, 935 (2009). 
i^D. Reinel and W. Dieterich, J. Chem. Phys. 104 5234 (1996). 
^^M. Kessler, W. Dieterich, H. L. Frisch, J. F. Gouyet, and P. Maass, Phys. Rev. E 65, 

066112 (2002). 
19 J. F. Gouyet, M. Plapp, W. Dieterich and P. Maass, Advances in Physics 52, 523 (2003). 
^°In the continuum hmit, and under the assumption of symmetric hopping rates this method 

has been shown^ to become equivalent with dynamic density functional theory^^^^ used 

in soft matter physics, see also^. 
21 W. Dieterich, H. L. Frisch and A. Majhofer, Z. Phys. B 78, 317 (1990). 
22U. Marini Bettolo Marconi and P. Tarazona, J. Chem. Phys. 110 8032 (1999). 
23p. Espaiiol and H. L5wen, J. Chem. Phys. 131, 244101 (2009). 
24A. B. Kolomeisky, Phys. Rev. Lett. 98, 048105 (2007). 

25A. M. Berezhkovskii and S. M. Bezrukov, Phys. Rev. Lett. 100, 038104 (2008). 
26A. Zilman, J. Pearson, and G. Bel, Phys. Rev. Lett. 103, 128103 (2009). 
2^M. Einax, M. Korner, P. Maass and A. Nitzan, Phys. Chem. Chem. Phys. 12, 645 (2010). 
2*M. Dierl, P. Maass, and M. Einax, in preparation. 
29b. Derrida, Phys. Rep., 1998, 301, 65 (1998). 
■^^G. Schiitz, in Phase Transitions in Critical Phenomena, edited by C. Domb and J. Lebowitz 

(Academic Press, San Diego, 2001), Vol. 19, pp. 3-251 
3iK. Sumithra and T. Sintes Physica A 297, 1 (2001). 
32j. K. Percus, Ace. Chem. Res. 27, 8 (1994) 

33 J. Buschle, P. Maass and W. Dieterich, J. Physics A 33, L41 (2000). 
3^If Pi + p2 > 1, a. situation that can be reached if fii and/or /i/j approach infinity at finite 

interaction V, a different solution is obtained as \^ — )■ oo that is not considered here. 
3^M. Einax, in preparation 

24 



36r. Evans, Adv. Phys. 28, 143 (1979). 

^^Since the numerator and denominator in Eq. fl36l) become zero for C = 0, the TDFT results 

labeled V = actually are for V = 0.001. 
38H. Spohn, J. Phys. A: Math. Gen. 16, 4275 (1983). 
39b. Derrida, J. Stat. Mech. P07023 (2007). 



25 






0.0125 


' (b) 






0.01 


- 


1 




0.0075 


- 


y 


1; 


0.005 
0.0025 


y^'^ 


% 









\ 






••A... 




-0.0025 


- 






-0.005 


- 






0.25 0.5 0.75 1 1.25 1.5 1.75 2 

t 



FIG. 1. Properties of a dilute system with pz, = p^ = 0.0067 (/-iL,R = —5.0), driving period r = 2 
and amplitude A = 5. The behavior of this dilute system does not depend on V for y = 0, 1 
and 10, and results from the TDFT approximation are indistinguishable from the exact ones, (a) 
Efficiency rj as a. function of mechanical bias F = er — e^ {ei = 0). Lower curve: neutral static site 
energies, e^ 2 = 0- Upper curve: attractive static site energies ^i 2 = —2.0. (b) Time dependent 
occupation probabilities pi(t), P2{t) and currents JL,R{t) for £12" —2.0 and F = 1.5. Shown is 
one period of stationary oscillation, where t = corresponds to a minimum in ei(i), see Eq. (jlip . 
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FIG. 2. (Color online) F-dependent efficiency r/ (a) and averaged current Jav (b) for different 
interparticle interactions 1/ in a driven channel coupled to moderately dense reservoirs with p/, = 
PR = 0.5 (/iL,_R = 0). Other parameters are e^ 2 = —2.0, r = 2, and A = 5.0. The insets show 
the V = Q case where interparticle interactions arise from site exclusion only, and the main figure 
shows results for V = —5, —1, 1, 10. Shown are results from the exact kinetic equations (red lines) 
and from the TDFT (blue lines). 
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FIG. 3. (Color online) Exact time dependent occupation probabilities Pi,2(i) (a)) left and right 
currents (b), as well as correlation function Pii{t) — Pi{t)p2it) for different y > (c and d) in a 
non-dilute system with pl,r = 0.5 and F = 1.5. Other parameters are as in Fig. [2l The plots 
(c) and (d) illustrate the differences between exact correlations and local equilibrium correlations 
assumed in the TDFT. 
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FIG. 4. (Color online) Efficiency rj (a) and averaged current Jav (b) from exact rate equations as 
a function of r for pl,b. = 0.5, V = 1 and different F = 0,1, 2, 3. Other parameters as in Fig. [2j 
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FIG. 5. (Color online) (a) The average current J^v from the exact rate equations in a peristaltic 
pump as a function of the electrochemical potential difference Afj,^ = F + A/i between the two 
reservoirs, for V = 1 and 10. Other parameters are as in Fig. El ^i 2 = —2.0, r = 2, and A = 5.0. 
Compared are the two limiting cases Jsx{F,0) where A/i = and Jav(0, A/z) where F = 0. In 
either case the left density is fixed to pL = 0.1 (/i^ « —2.197). (b) Same for a flashing ratchet with 
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-1.0, r = 2, and A = 2.0. 
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FIG. 6. (Color online) Chemical efficiency {F = 0) compared with mechanical efficiency (A^ = 0) 
for peristaltic driving, both as a function of A/.i^ = F + A^ and for different V. The left density 
is fixed to pi = 0.1 {m ~ —2.197). Other parameters as in Fig. [5^. 
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FIG. 7. (Color online) Efficiency (a) and average current (b), plotted against the bias F for a 
flashing ratchet in the moderately dense system characterized by pL = Pr = 0.5 {^l.r = 0). Other 
parameter are £-[ 2 = —1-0, r = 2, and A = 2.0. 
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FIG. 8. (Color online) The ideal efficiency rjs displayed as a function of the interparticle coupling 
V using T = 2 (a) and of the modulation time r using V = 1 (b), for a peristaltic pump with 
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-2.0, A = 5 and three different choices for calculating Win- at F = 0, J = and at 



F = i^max) where the efficiency rj defined by Eq. ([9]) is maximal. 



33 



